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Abstract 

We investigate spherical macroions in the strong Coulomb coupling regime within 
the primitive model in salt-free environment. Molecular dynamics (MD) simulations 
are used to elucidate the effect of discrete macroion charge distribution on charge 
inversion. A systematic comparison is made with the charge inversion obtained in the 
conventional continuous charge distribution. Furthermore the effect of multivalent 
counterions is reported. 
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1 Introduction 



Charged colloids (macroions) might exhibit very surprising behaviors due to 
the Coulomb interaction and the two different length scales involved there: (i) 
microscopic (counterions) and (ii) mesoscopic (macroions) . A very spectacular 
effect which has attracted great attention these last years is the phenomenon 
of overcharging [1-8]. This situation occurs when the number of counterions 
in the vicinity of the macroion surface is so high that the macro-particle bare 
charge is overcompensated (charge inversion). This phenomenon has been ob- 
served experimentally by electrophoresis [9,10]. On the theoretical side, a mas- 
sive effort has been recently devoted to elucidate this striking phenomenon. 
Wigner crystal theories have been successfully applied to understand charge 
inversion [1,3-6,8] and quantitative agreement was found with simulations [8]. 

It turns out that all numerical and analytical methods, undertaken so far, 
neglect the discrete nature of the macroion charge distribution. But in re- 
alistic systems the macroion charge is not continuously distributed over the 
macro-particle, but rather carried by small microscopic ions. It is only very 
recently that a simulation study has been reported where the macroion charge 
discreteness has been considered [7]. 
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In this paper, we concentrate on the effect of the macroion charge distribution 
discreteness. In particular, we study this effect on the overcharging properties. 
Besides we show that it is possible to have stronger overcharge in the discrete 
case than in the continuous case. 



2 Computational details 

2.1 Macroion charge discretization 

The procedure is similar to the one used in a previous study [7] . The macroion 
charge discretization is produced by using Z rn monovalent microions of diam- 
eter a distributed randomly on the surface of the macroion (see Fig. 1). Then 
the macroion structural charge is Q = —Z m e [11] where Z m > and e is the 
positive elementary charge. The counterions have a charge q = +Z c e where 
Z c > stands for the counterion valence. The discrete colloidal charges (DCC) 
are fixed on the surface of the spherical macroion and are at a distance r from 
the macroion center. Thus to produce a random discrete charge distribution 
on the surface we generated (uniformly) randomly the variables cos 9 and <p 
where 9 and ip are the standard spherical angle coordinates. Excluded vol- 
ume is taken into account by rejecting configurations leading to an overlap 
of microions. Phenomena such as surface chemical reactions [12], hydration, 
roughness [13] are not considered. 



2.2 Molecular dynamics procedure 

A MD simulation technique was used to compute the motion of the counterions 
coupled to a heat bath acting through a weak stochastic force W(t). The 
motion of counterion i obeys 

cP"v- dv ■ 

m—^ = -Vitffr) - m 7 -^ + W 4 (f) , (1) 



where m is the counterion mass, U is the potential force having two contri- 
butions: (i) the Coulomb interaction and (ii) the excluded volume interaction, 
and 7 is the friction coefficient. Friction and stochastic force are linked by the 
dissipation-fluctuation theorem < Wj(i) • ~Wj(f) >= 6m , ykBTSij5(t — t). For 
the ground state simulations the stochastic force is set to zero. 
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Excluded volume interactions are taken into account with a pure repulsive 
Lennard- Jones potential given by 
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(2) 



if r — r < r cut , and Ulj(t) = otherwise. For the microion-microion interac- 
tion (the microion being either a counterion and/or a DCC), whereas ro = 7cr 
for the macroion-counterion interaction, and r cut = 2 1 / 6 cr is the cutoff radius. 
This leads to a macroion-counterion distance of closest approach a = 8cr. 

Energy and length units in our simulations are related to experimental units 
by taking e =A;bT (with T = 298 K), and a = 3.57 A respectively. 

The pair electrostatic interaction between any pair ij, where i and j denote 
either a DCC, a counterion or the central charge (for the case of the continuous 
surface charge distribution), reads 



where Ib = e 2 /4neoe r kBTo is the Bjerrum length (here Ib = 10a) describing 
the electrostatic strength. 

The macroion and the counterions are confined in a spherical impenetrable cell. 
The macroion is held fixed and is located at the center of the cell. To avoid 
image charge complications, the permittivity e r is supposed to be identical 
within the whole cell (including the macroion) as well as outside the cell. 



3 Results 

The simplest way to quantify overcharging in salt-free environment is to com- 
pute the total electrostatic energy of the system in the ground state (i. e., 
T = OK) as a function of the number of overcharging counterions n. Although 
this approach is very idealized, it allows an insight into the mechanisms gov- 
erning overcharging. 

In what follows we are going to compare systematically (for a macroion charge 
Z m = 60) the overcharging occurring in the continuous case (CC), as is conven- 
tionally done, and in the discrete case (with DCC ions). Fig. 1 shows a typical 
starting equilibrium configuration with a discrete macroion charge distribu- 
tion where the system is globally neutral (before adding excess overcharging 
counterions) . 




(3) 
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Fig. 1. Ground state structure for a discrete macroion charge distribution with 
Z m = 60. The monovalent discrete colloidal charges are in white and the monovalent 
counterions (here Z c = 1) in blue. Full ionic pairing occurs. 

Fig. 2(a) shows overcharging energy profiles for Z m = 60 and monovalent 
counterions (Z c = 1). It is observed that overcharging is stronger in the con- 
tinuous case than in the discrete case as one could expect. More precisely 
the gain in energy as well as the maximal number of stabilizing overcharging 
counterions (corresponding to a minimum in the E — n curve) are higher when 
the surface macroion charge is continuous. Indeed, having in mind that over- 
charging is enhanced by counterion ordering, it is clear that in the discrete 
case where strong ionic pairing occurs (between DCC site and counterion - see 
Fig. 1) the counterion ordering is lowered and thus leads to a weaker charge 
inversion. However since overcharging occurs in the discrete case, the excess 
counterions still experience an attractive interaction with their neighboring 
dipoles (ionic pairs). This ion-dipole interaction is the driving force for the 
overcharging in the discrete case with monovalent counterions. Nevertheless 
it has been shown elsewhere [7] that by decreasing the mean separation be- 
tween DCC sites (i. e., increasing macroion charge density) the difference in 
overcharging between discrete and continuous macroion charge distribution 
decreases. 

When we consider multivalent counterions the situation is qualitatively differ- 
ent. Indeed, one observes that for Z c = 2 and 3 overcharging is even stronger 
in the discrete case [see Figs. 2 (b-c)]. This is a rather counter-intuitive phe- 
nomenon, since ionic pairing is still effective. The crucial difference between 
monovalent and multivalent counterion systems is that overcharging occurs 
with ionic pairing for multivalent counterions (since free monovalent DCC sites 
are unbound), whereas for monovalent counterions only ion-dipole interaction 
can take place (since in the neutral state each DCC site is already bound with 
a monovalent counterion). Obviously for these systems under consideration 
depicted in Figs. 2 (b-c), the ionic pairing energy is a fundamental quantity 
that controls the overcharging strength. More precisely, the stronger the ionic 
bonding the more efficiently it overcomes the excess overcharge " self-energy" , 
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Fig. 2. Total electrostatic energy for ground state configurations (with Z m = 60) as 
a function of the number of overcharging counterions n: (a) Z c = 1 (b) Z c = 2 (c) 
Z c = 3. Overcharge curves were computed for discrete macroion charge distribution 
(DCC) and central macroionic charge (CC). The neutral case was chosen as the 
potential energy origin. Dashed lines were produced by using Wigner Crystal theory 
[8]. For discrete systems (DCC) error bars are only indicated when larger than 
symbols. 



and therefore the higher the overcharging for multivalent counterion systems. 

Beside of the ionic bonding strength there is another important parameter con- 
trolling overcharging: the surface macroion charge density. When the macroion 
charge density is high enough (not reported here) then the E(n) curves be- 
come quasi-identical for discrete and continuous charge distributions. In terms 
of length scales, the relevant parameter controlling charge reversal with dis- 
crete distribution is the ration p between the mean distance between DCC and 
the distance separation of the ionic pair. For low p the effects of discretization 
are very important whereas for p — > 1 one recovers the continuous limit. 
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4 Conclusions 



We carried out MD simulations to elucidate the effect of discrete macroion 
charge distribution on the overcharging. We showed that overcharging is still 
effective in the discrete case. For multivalent counterions and a sufficiently low 
macroion charge density, overcharging can even be stronger in the discrete case 
than in the continuous case. Future works will address aqueous solutions at 
finite temperature as well as the presence of salt-ions. 
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